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Abstract 


A search is described for a Higgs boson decaying into two photons, one of which 
has an internal conversion to a muon or an electron pair (M'y). The analysis is per¬ 
formed using proton-proton collision data recorded with the CMS detector at the 
LHC at a centre-of-mass energy of 8 TeV, corresponding to an integrated luminos¬ 
ity of 19.7fb^^. The events selected have an opposite-sign muon or electron pair and 
a high transverse momentum photon. No excess above background has been found 
in the three-body invariant mass range 120 < niwy < 150 GeV, and limits have been 
derived for the Higgs boson production cross section times branching fraction for the 
decay H —>• 7*7 —)■ iij, where the dilepton invariant mass is less than 20 GeV. For 
a Higgs boson with wth = 125 GeV, a 95% confidence level (GL) exclusion observed 
(expected) limit is 6.7 (5.9l^ g) times the standard model prediction. Additionally, an 
upper limit at 95% GL on the branching fraction of H —)■ for the 125 GeV Higgs 

boson is set at 1.5 x 10^^. 
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1 Introduction 

The rare decay into the ii'y final state of the Higgs boson is a rich source of information that 
can enhance our understanding of its basic properties and probe novel couplings predicted by 
extensions of the standard model (SM) of particle physics. As illustrated in Fig. this decay 
in SM has contributions from loop-induced H —> 7*7 and H —> Z7 diagrams (a, b, c), tree- 
level process H —)> with final-state radiation (d), and higher-order processes, known as box 

diagrams (e, f, g) UMI. Other contributions include H —> V{qq)j —)■ shown in Fig. 

where V denotes a vector meson Q/xp or Y) that decays to £i BQ- The Higgs boson branching 
fraction to ££'y is dominated by the H —> 7*7 and H —> Z7 modes, while the contribution from 
the box diagrams is negligible [ 1 ]. In the muon channel, when the dilepton invariant mass, 
mi£, is greater than 100 GeV, final-state radiation in H ^ fifi starts to dominate [[ 81 . In the 
three-body decay, H ^ ££j, it is possible to investigate non-SM couplings by examining the 
angular distributions, and forward-backward asymmetry variables reconstructed from the ££'y 
final state HU 13 . 

The expected rates of the H —> {Z/j*)'y —> ££'y processes compared to the rate of H ^ 77 
decay, for a Higgs boson with mass mu = 125 GeV, are lUOUTTI : 


r(H —> 7*7 —>• 667) 
r(H^77) 


3 . 5 %, 


r(H —>■ 7*7 —> jipiy) 
r(H ^ 77) 


1.7%, 


r(H ^ Z7 ^ Hj) 

r(H^77) ^ 


2.3%. 



Figure 1 : Diagrams contributing to H —)■ ££j. The contributions from diagrams (a), (b), and 
(c) dominate. The final-state radiation of H —> decay, shown in diagram (d), is important 

at high dilepton invariant mass. Higher order contributions from diagrams (e), (f) and (g) are 
negligible. 


The H —)■ 7*7 —)■ ee7 decay is distinct from H —> 77 followed by a conversion of a photon to an 
e+e^ pair in the detector, which can become a background for H —)■ 7*7 if photon conversions 
are not properly identified. Experimentally, the various contributions shown in Figs. [^andll 
can be disentangled to some extent. Requirements on and the transverse momentum (px) 
of the photon are used to separate H —> 7*7 and H —> Z7. Events with final-state radiation are 
removed by requiring the photon to be isolated from either of the leptons. Gontributions from 
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2 CMS detector and trigger 



H ^ ^ ^^7 and other resonances are identified and rejected or selected based on the 

value of m££. 

The ATLAS and CMS Collaborations at the CERN LHC have both performed a search for 
H ^ Z 7 —> ii'Y decay with above 50 GeV [l2j[T3l. As a natural extension of those analyses, 
the current paper describes the first search for a Higgs boson Dalitz decay, H —> 7 * 7 , where 
the 7 * decays into a muon or an electron pair. The search is performed for a Higgs-like particle 
within the mass range between 120 and 150 GeV. In order to select the contribution from the 
Dalitz decay, we require niu < 20 GeV. The topology is a clean final state with a mass 
resolution of about 1.6%, as measured from the simulated signal samples. The eey channel is 
challenging due to the low ni££ that results in a pair of merged electron showers in the electro¬ 
magnetic calorimeter (EGAL). Nevertheless, when the merged showers are reconstructed in the 
EGAL, a mass resolution of 1.8% is achieved. Important backgrounds include the irreducible 
contributions from the initial- and final-state photon radiation in Drell-Yan production, and 
Drell-Yan events with additional jets where a jet is misidentified as a photon. 

In addition, a search is performed for H —)■ —?■ jiji'y decay for niy^ = 125 GeV, which 

is sensitive to the Higgs boson coupling to charm quark and a promising way to access the 
couplings of the Higgs boson to the second generation quarks at the LHG. In the SM this decay 
occurs through two main processes: direct coupling of the Higgs boson to charm (Eig.[^), and 
the usual t/W loop, where the radiated 7 * converts to a cc in a resonant state (Eig. |^). The 
two amplitudes interfere destructively and the second one dominates lUlZl- Eor the SM Higgs 
boson with wjh = 125 GeV, the branching fraction is predicted to be 2.8 x 10^^. A search by the 
ATLAS Gollaboration for this decay is described in Ref. IMH . 

The results presented in this paper are based on proton-proton collision data recorded in 2012 
with the GMS detector at a centre-of-mass energy = 8 TeV, corresponding to an integrated 
luminosity of 19.7 fb^^. 

2 CMS detector and trigger 

A detailed description of the GMS detector, together with a definition of the coordinate system 
used and the relevant kinematic variables, can be found in Ref. Ill5l . The central feature of 
the GMS apparatus is a superconducting solenoid, 13 m in length and 6 m in diameter, which 
provides an axial magnetic field of 3.8 T. Within the solenoid volume are a silicon pixel and 
strip tracker, the EGAL, and a hadron calorimeter (HGAL). Gharged-particle trajectories are 
measured by silicon pixel and strip trackers, covering 0 < (p < 2n in azimuth and \ri\ < 2.5 in 
pseudorapidity. A lead tungstate crystal EGAL surrounds the tracking volume. It is comprised 
of a barrel region \rj\ < 1.48 and two endcaps that extend up to |?/1 = 3. A brass and scintillator 
HGAL surrounds EGAL and also covers the region \r]\ < 3. Iron forward calorimeters with 
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quartz fibers, read out by photomultipliers, extend the calorimetric coverage up to \ri\ = 5. 
A lead and silicon-strip preshower detector is located in front of the ECAL endcaps. Muons 
are identified and measured in gas-ionization detectors embedded in the steel flux-return yoke 
outside the solenoid. The detector is nearly hermetic, allowing energy balance measurements 
in the plane transverse to the beam direction. 

A two-tier trigger system selects collision events of interest for physics analysis. Two triggers 
are used in the current analysis. In the muon channel, the trigger requires a single muon and 
a photon, both with pj greater than 22 GeV. In the electron channel the 7 * —> ee process at 
low dielectron invariant mass mimics a photon at the trigger level. For this reason, a diphoton 
trigger is used in the electron channel, for 7 -I 7 * final state. The trigger requires a leading 
(subleading) photon with pj > 26 (18) GeV. The diphoton trigger is inefficient for events with 
high dielectron invariant mass (nigg > 2 GeV) due to the isolation and shower shape require¬ 
ments. The available dielectron triggers cannot be used to select events with 2 < mgg < 20 GeV 
because they also require isolation, and the pj requirement made on the subleading lepton is 
too high. 

3 Event reconstruction 

The photon energy is reconstructed from a sum of signals in the EGAL crystals IIT^ . The EGAL 
signals are calibrated and corrected IflTl , and a multivariate regression technique, developed for 
the H —> 77 analysis Ill 8 l . is used to determine the final energy of the photon fl^ . The neigh¬ 
boring EGAL crystals with energy deposition are combined into clusters, and the collection of 
clusters that contain the energy of a photon or an electron is called a supercluster. Identification 
criteria are applied to distinguish photons from jets and electrons. The observables used in the 
photon identification criteria are: the isolation variables, the ratio of the energy in the HGAL 
towers behind the supercluster to the electromagnetic energy in the supercluster; the trans¬ 
verse width in p of the electromagnetic shower; and the number of charged tracks matched to 
the supercluster. The efficiency of the photon identification is measured using Z —> ee data by 
reconstructing the electron showers as photons, and found to be 80 ( 88 %) at a transverse energy 
> 30 (50) GeV and \p\ < 1.44. 

Muon candidates are reconstructed in the tracker and identified by the particle-flow global 
event reconstruction algorithm IIT^ using hits in the tracker and the muon systems. This 
approach allows us to maintain a high efficiency independent of the dimuon invariant mass 
and to reconstruct muons with pj as low as 4 GeV. Muons from ^ pp internal conversions 
are expected to be isolated from other particles. A cone of size AR = (A//)^ -|- (A^)^ = 0.4 is 

constructed around the momentum direction of each muon candidate 1211 . The relative isola¬ 
tion of the muon is quantified by summing the pj of all photons, charged and neutral hadrons 
within this cone, and then dividing by the muon pj. The resulting quantity, corrected for ad¬ 
ditional underlying event activity due to pileup events, is required to be less than 0.4 for the 
leading muon. The isolation requirement rejects misidentified leptons and background arising 
from hadronic jets. The AR{pp) separation between the two muons is small due to their low 
invariant mass (as shown in Fig.|^ and high pj of the 7 * in H —> 7*7 decays. Hence, no iso¬ 
lation requirement is applied to the subleading muons as they are already within the isolation 
cone of the leading muons in most events. Dimuon identification and isolation efficiency of 
about 80% is obtained. 

In the electron channel of the H —> 7*7 —)> decay, the two electrons produced in the 7 * —> ee 

process are even closer to each other than in the muon channel, since the niu is smaller (Fig.|^. 
Therefore, their energy deposits in the EGAL are merged into one supercluster giving rise to 
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5 Event selection 


a unique signature. To identify these merged electrons, two tracks associated to the super¬ 
cluster are required. A Gaussian sum filter (GSF) algorithm is used to reconstruct the electron 
tracks Il22l . The supercluster energy must correspond to pj > 30GeV and be located in the 
EGALbarrel {\r]\ < 1.44). The scalar sum pj^ + p^ of the corresponding two GSF tracks must 
exceed 44 GeV. Both GSF tracks are required to have no more than one missing hit in the pixel 
detector in order to reduce the background from photons converting to e+e^ in the detector 
material. A multivariate discriminator is trained to separate the 7 * —> ee objects from jets or 
single electrons. The input variables for the training include lateral shower shape variables, 
the median energy density in the event to take into account the pileup dependence, and the 
kinematic information from the supercluster and tracks. A combined reconstruction and selec¬ 
tion efficiency of ~40% is achieved for the signal. For comparison, the efficiency for a single 
isolated electron with similar pj is ~ 88 % Il23ll . 


4 Simulated samples 

The description of the Higgs boson signal used in the search is obtained from simulated events. 
The samples for the Dalitz signal are produced at leading-order using the MadGraph 5 
matrix-element generator Il24l with the ANO-HEFT model ||25l , interfaced with pythia 6.426 
|j26l , for the gluon and vector boson fusion processes, and for associated production with a 
vector boson. Associated production with a tf pair is ignored because of its small contribution. 
The sample for H —> process is produced with the pythia 8.153 generator I27l, and 

reweighted to simulate 100% polarization of the p/xp. The parton distribution function (PDF) 
set used to produce these samples is given by GTEQ6L1 Il28l . The SM Higgs boson production 
cross sections are taken from Ref. IITT] . The branching fractions for H —> 7*7 are estimated 
using MGFM 6.6 Il29ll and for H —> {]/ip)'y are taken from Ref. DU. For the SM Higgs boson in 
the mass range of 120-150GeV, the H —>• 7*7 —>• ppj (ee'y) branching fraction is expected to 
be between 2.0 (4.5) x 10^^ and 3.3 (7.5) x 10^^ for below 20 GeV. The expected branching 
fraction for H —>• {]/ip)'Y is (2.8 ± 0.2) x 10^^ for ?Wh = 125 GeV, which is further suppressed 
due to the ]/xp meson decay to muons, BQ/xp —> pp) = 0.059. 

The simulation aims to include all known effects and the conditions of real data taking in GMS. 
Some residual differences between the data and simulation are taken into account by reweight¬ 
ing the simulated events with scale factors. Systematic uncertainties are assigned to cover 
imperfect knowledge of residual differences. Scale factors are implemented to match the dis¬ 
tribution of primary vertices, the photon identification and isolation efficiency, and the muon 
isolation efficiency. No corrections are applied to the muon and electron identification and 
trigger efficiencies, but an uncertainty is assigned as described in Section]^ 

The energy and momentum resolution of muons and photons in simulated events are corrected 
to match that in data. The energy scale of muons (photons) is corrected to that found in Z —> 
pp (ee) events. For the electrons, no resolution or scale corrections are applied because of 
their unique topology, and the absence of a data-driven method to derive those corrections. 
Therefore, we rely on the simulation of the 7 * —> ee process and assign uncertainties sufficient 
to cover any possible discrepancy in the scale and resolution between data and simulation. 


5 Event selection 

Events are required to pass the muon plus photon trigger in the ppj final state and the dipho¬ 
ton triggers in the eej final state. The trigger efficiency for signal events after the selection 
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Figure 3: The invariant mass of the dilepton system in signal simulation for nin = 125 GeV. 
Distributions are shown for muon and electron channels, before and after selection. The in¬ 
variant mass before selection is obtained from the leptons at the generator level, while after 
selection the reconstructed invariant mass is used. 

requirements described below is 85% (90%) in the muon (electron) channel, as measured from 
the simulated samples. 

The muons (electrons) are required to be within \ri\ < 2.4 (1.44), while the photon is required 
to be within \ri\ < 1.44. The invariant mass of the ££7 system, is required to satisfy 

110 < m££y < 170 GeV. The photon and dilepton momenta both must satisfy pj > 0.3 ■ m££y 
requirement, which is optimized for high signal efficiency and background rejection. 

On average, there are 21 pp interactions within the same bunch crossing in the 8 TeV data, 
which result in about 16 collision vertices reconstructed in each event. The vertex with the 
highest scalar sum of the pj of its associated tracks is taken to correspond to the primary inter¬ 
action vertex. The primary vertex must have the reconstructed longitudinal position (z) within 
24 cm of the geometric centre of the detector and the transverse position (x-y) within 2 cm of the 
beam interaction region. The lepton tracks from 7 * —> pp (ee) are required to originate from 
the primary vertex, and to have transverse and longitudinal impact parameters with respect to 
that vertex smaller than 2.0 (0.2) mm and 5 (1) mm, respectively. 

The muons must be oppositely charged, and have pj > 23 (4) GeV for the leading (subleading) 
lepton. The pj requirement on the leading muon is driven by the trigger threshold, and on the 
subleading muon by the minimum energy needed to reach the muon system, while maintain¬ 
ing high reconstruction efficiency. In the electron channel, no additional selection on pj of the 
GSF tracks is necessary, beyond those described in Section]^ Finally, in both muon and electron 
channels, the separation between each lepton and the photon is required to satisfy AR > 1 in 
order to suppress Drell-Yan background events with final-state radiation. 

The dilepton invariant mass in the muon channel is required to be less than 20 GeV to reject con¬ 
tributions from pp —> 7 Z and to suppress interference effects from the H —> 7 Z process and the 
box diagrams shown in Fig. Events with a dimuon mass in the ranges 2.9 < < 3.3 GeV 

and 9.3 < < 9.7 GeV are rejected to avoid the ]/ip —> pp and Y ^ pp contamination. In 

the electron channel the invariant mass, constructed from the two GSF tracks, is required to 
satisfy wJee < 1.5 GeV. The m££ distributions for simulated signal events are shown in Fig. |^in 
the muon and electron channels. 
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6 Background and signal modeling 




Figure 4: The (left) and nteej (right) spectra for 8 TeV data (points with error bars), together 
with the result of a background-only fit to the data. The Icr and 2cr uncertainty bands represent 
the uncertainty in the parameters of the fitted function. The expected contribution from the SM 
Higgs boson signal with mn = 125 GeV, scaled up by a factor of 10, is shown as a histogram. 


In the search for the H —> both pj and pj^ > 40GeV are required, and the 

events are selected with 2.9 < < 3.3 GeV. 

The observed yields after the event selection described above are listed in Table In the elec¬ 
tron channel, there is also a contribution from the H —> 77 process due to unidentified conver¬ 
sions, which is about 15% of the H —> 7*7 signal (0.2 events at wjh = 125 GeV). This contri¬ 
bution is considered as a background to H —)■ 7 * 7 , and negligible compared to the continuum 
background estimated from the fit to data described in the next section. 

Table 1: The expected signal yield and the number of events in data, for an integrated luminos¬ 
ity of 19.7fb^^. Signal events are presented before and after applying the full selection criteria 
described in the text. In the (J/t/’)7 sub-category only the ]/tp pp decay is considered, and 
the signal yield is a sum of two contributions: H —> > ppj and H —> 7*7 —> ppj, 

where the dimuon mass distribution is non-resonant. 



Signal events 

Signal events 

Number of events 

Sample 

before selection 

after selection 

in data 


oth = 125 GeV 

mn = 125 GeV 

120 < m££y < 130 GeV 

ppX 

13.9 

3.3 

151 

eej 

25.8 

1.9 

65 

Q/ip pp)x 

0.065 (J/i/ 7 ) -|- 0.32 (non-res.) 

0.014(J/i/;) - 1 - 0.078 (non-res.) 
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6 Background and signal modeling 

The background is modeled by fitting a polynomial function to the ££'y mass distributions in 
data. An unbirmed maximum likelihood fit is performed over the range 110 < ni££j < 170 GeV. 
Fig.|^shows the m££y spectra, which are fitted with pol}momial functions of fourth degree. The 
reduced of tho fits are 0.5 and 0.7 for the muon and electron channels, respectively. Even 
though the search is limited to 120 < wjh < 150 GeV, the fits to the m££y spectra are performed 
over a wider range, giving a better modeling of the background, particularly at the edges of the 
search range. The degree of the pol}momials is chosen following a procedure similar to the one 
described in Ref. I30l . This procedure ensures that the potential bias due to the background 
modeling is at least five times smaller than statistical uncertainty. 

For the H —> search, where only the single Higgs boson mass hypothesis jhh = 125 GeV 























7 


is investigated, a fit to a polynomial of second degree is performed over the 110-150 GeV mass 
range (Fig.|^. 

The signal model in all three cases is obtained from an unbinned fit to the mass distribution of 
the corresponding sample of simulated events to a Crystal Ball function plus a Gaussian 
function. 
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Figure 5: The distribution for events with 2.9 < < 3.3 GeV for 8 TeV data (points 

with error bars), together with the result of a background-only fit to the data. The Icr and 
2cr uncertainty bands represent the uncertainty in the parameters of the fitted function. The 
expected contribution from the H —> (J/ip)j —> process of the SM H with = 125 GeV, 

scaled up by a factor of 500, is shown as a histogram. 


7 Results 

The data are used to derive upper limits on the Higgs boson cross section times branching 
fraction, C 7 (pp ^ H) S(H —> 7*7 —> ££j) divided by that expected for a SM Higgs boson, 
for m££ < 20 GeV. No significant excess above background is observed in the full mass range, 
120 < ?Wh < 150 GeV, with a maximum excess of less than two standard deviations. In the elec¬ 
tron channel a correction is made to account for the events that are removed by the requirement 
of nigg < 1.5 GeV due to the trigger and reconstruction inefficiencies described above. 

The exclusion limits are calculated using the modified frequentist GLg method Il32] - l36l . An 
unbrnned evaluation over the full mass range of data is used. The uncertainty in the limit is 
dominated by the size of the data sample and systematic uncertainties have a small impact. 

The systematic uncertainty in the limits results only from the uncertainty in the signal descrip¬ 
tion, as the background is obtained from data and biases in the fitting procedure have been 
found to be negligible. A summary of the systematic uncertainties is given in Table The un¬ 
certainty can be separated into the uncertainty resulting from theoretical predictions and from 
the uncertainty in detector reconstruction and selection efficiency. 

Theoretical uncertainties come from the effects of the PDF choice on signal cross section, the 
missing higher-order calculations (scale) Il38ll42ll . and the uncertainty in the prediction on the 
Higgs boson decay branching fraction (Ullll. The uncertainty due to the muon reconstruction 
efficiency, 11%, is obtained from data using ]/tp —> }i}i events. It is dominated by the statistical 
uncertainty of the data sample. In the electron channel, the corresponding uncertainty, 3.5%, is 
obtained from simulation. The 11% uncertainty estimated for the muon identification efficiency 
is sufficiently small and it has no impact on our result, thus no simulation study was attempted, 
although it could greatly reduce the uncertainty. 
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8 Summary 


Table 2: Systematic uncertainties affecting the signal 


Source 

Uncertainty 

Integrated luminosity (ref. II37I) 

2 . 6 % 

Theoretical uncertainties: 


PDF 

2.6-7.5% 

Scale 

0.2-7.9% 

H —> 7*7 —>• ££7 branching fraction 

10 % 

Experimental uncertainties: 


Pileup reweighting 

0 . 8 % 

Trigger efficiency, /r (e) channel 

4(2)% 

Muon reconstruction efficiency 

11 % 

Electron reconstruction efficiency 

3.5% 

Photon reconstruction efficiency 

0 . 6 % 

m£i,y scale, }i (e) channel 

0.1 (0.5)% 

m££y resolution, (e) channel 

10 ( 10 )% 


The expected and observed individual and combined jifi'y and eey limits are shown in Fig. 
The limits are calculated at 1 GeV intervals in the 120-150 GeV mass range. The median ex¬ 
pected exclusion limits at 95% confidence level (GL) are between 6 and 10 times the SM pre¬ 
diction and the observed limit ranges between about 5 and 11 times the SM. The observed 
(expected) limit for niu = 125 GeV is 6.7 (5.9^^ g) times the SM prediction. 

The 95% GL exclusion limits on cr(pp —)■ H) 13(H —> for a narrow scalar particle without 

assuming the decay kinematics of a SM Higgs boson, in the muon channel, are shown in Fig. 
The observed (expected) limit for wjh = 125 GeV is 7.3 (5.2^^ ^) fb. The total signal efficiency is 
24% and almost independent of the dimuon invariant mass. In the electron channel, however, 
this efficiency depends on the dielectron mass, since it is strongly shaped by the selection. For 
this reason the corresponding limit in the electron channel is not evaluated. 

Additionally, for the SM Higgs boson with ?mh = 125 GeV, we place an upper limit for a 2.9 < 
ntii < 3.3GeV region in the muon channel: cr(pp —> H)13(H —> < 1.80 fb, while the 

expected limit is 1.90 ± 0.97 fb. One can interpret this result as an upper limit on cr(pp —> 
H) 13(H ^ (J/V’)7 FFT) ^rid obtain for the branching fraction, 5(H —> < 1.5 x 10^^ 

at 95% GL, which is about 540 times the prediction in Ref. ||i6i|. The limit on the branching 
fraction at 90% GL is 5(H —> < 1.2 x 10“^. The number of events present in this 

mass window coming from the H —)■ 7*7 —)■ /r ^7 is large compared to the H —)■ (J/j /’)7 —>• 

(as shown in Table [^. On the other hand it is small compared to the total background, hence it 
is considered as a part of the background when extracting the limit on 6 (H —> 


8 Summary 

A search for a Higgs boson decay H —> 7*7 —)■ ££^ is presented. No excess above the back¬ 
ground predictions has been found in the three-body invariant mass range 120 < m££y < 
150 GeV. Limits on the Higgs boson production cross section times the H —)■ 7*7 — 
branching fraction divided by the SM values have been derived. The observed limit for wjh = 
125 GeV is about 6.7 times the SM prediction. Limits at 95% GL on cr(pp ^ H) S(H —> 
for a narrow resonance are also obtained in the muon channel. The observed limit for wjh = 
125 GeV is 7.3 fb. Events consistent with the ]/tp in dimuon invariant mass are used to set a 
95% GL limit on the branching fraction 13(H —> (J/t/’) 7 ) < 1.5 x 10^^, that is, 540 times the SM 
prediction for mu = 125 GeV. 
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Figure 6: The 95% CL exclusion limit, as a function of the mass hypothesis, mu, on a/ ctsm/ the 
cross section times the branching fraction of a Higgs boson decaying into a photon and a lepton 
pair with m££ < 20 GeV, divided by the SM value. (Upper left) muon, (upper right) electron 
channels, (bottom) statistical combination of the results in the two channels 
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8 Summary 



Figure 7: The 95% CL exclusion limit on (7(pp —> H) i3(H ^ FFT)/ with < 20GeV, for a 
Higgs-like particle, as a function of the mass hypothesis, ?«h- 
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